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INTRODUCTION 


A single component laser Doppler velocimeter (LDV) system provides a 
technique for measuring one dimensional mean and turbulent velocities in 
fluid flows. It is shown here, that the single component LDV system also 
is capable of determining the additional information of two and three dimen- 
sional mean and turbulent flow characteristics. This capability is illustra- 
ted by the derivation of several sets of equations which define the mean and 
turbulent parameters that a single component LDV system measures in a three 
dimensional flow field. It is shown that by making LDV measurements from sev- 
eral different independent angular arrangements, sets of independent simultaneous 
equations are defined and can be solved to obtain three dimensional mean, RMS, 
cross-products and correlation parameters of the fluid flow. This procedure 
assumes that the independent measurements are made at the same location in the 
flow and that the flow is statistically stationary. 

Solutions to the sets of simultaneous equations were considered for a 
single component forward scatter LDV system, to assess the precision with which 
the Doppler frequencies must be measured in order to obtain an accurate measure- 
ment of the three dimensional flow paramenters. It was found that for the ang- 
ular arrangements considered, the three dimensional mean flow quantities would 
be the least difficult measurement, while the RMS, cross-products and correla- 
tion parameters become more difficult. This difficulty is due to a more string- 
ent requirement imposed upon the precision to which the Doppler frequencies 
must be evaluated. 



LDV SYSTEM AND DEFINING EQUATIONS 

The laser Doppler velocimeter employs the phenomena of Doppler shifting 
of light frequencies to measure the velocities of small particles suspended 
in fluid flows. For nonrelativistic particle motion, the Doppler shift may be 
shown to be linearly related to the velocity of a particle in a direction 
dependent on the geometrical angular arrangement of the optical paths within 
the system. Figure 1 shows three typical optical arrangements that are used 
in single LDV systems. For each of these systems a general equation can be 
written which relates the instanstaneous Doppler frequency measured by the 
system to the instanstaneous velocity components of a cartesian coordinate 
system. This equation is written as 

XF.(t) = a. U(t) + b. V(t) + c. W(t) (d 

Where U, V, and W are the velocity components in the X, Y, and Z directions, 

A is the laser wavelength, F^ is the Doppler frequency shift, a^ , b ^ , c.. 
are coefficients relating the geometrical arrangement of the LDV system to a 
cartesian coordinate system, and the subscript i denotes the measurement 
number for a given system configuration, i.e. i = 1,2, n. 

Dividing each of the velocities and the Doppler frequency into a mean 
and a fluctuating component, such as: 

U(t) = UftJ +■ U'(t) (2) 

Equation 1 becomes: 

X(Mt) + F, 1 (t ) ) = aJUft) + U'(t)l + b [v(tj + V'(t)J 

1 c, [«t> ♦ WU)] ’ (3) 

The overbar indicates the quantities mean value or time average, which may 

be expressed mathematically as: 
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(4) 


( ) = _L / ( )dt 

T o 

and from the definition in (2) ()'(t) 0 

Taking the time average of Equation 3 yields: 

X F. (t) = a. U ( t j + b 1 T{t) + Ci W (t) (5) 

Where U(t), V (t ) , and W(t) are the mean velocities in the X, Y, and Z directions 
and F..(t) is the mean Doppler frequency shift for a single LDV system in a 
configuration defined by the a^ , , and c^ coefficients. 

From Eq. 5 it is seen that in a general single LDV measurement, there 
are three unknown velocity components. These velocities may be determined by 
making LDV measurements at the same position in the flow field and changing 
the geometrical arrangement of the system for three independent measurements. The 
defining equations for such a set of measurements are written as: 

X F 1 (t) = a, uiT) + b l ~Ht) + CjIkT) (6) 

X F 2 (t) = a 2 ~U(tj + b 2 "v(t) + c 2 W(t) (7) 

X F 3 (t) = a 3 U(tj + b 3 "v(tj + c 3 W(t) (8) 

By solving Eqs. (6), (7), and (8) simultaneously the three mean velocity 
components may be obtained. 

Since the three indicated measurements must be conducted at three different 
times, the flow must be stationary in order for the procedure to yield meaning- 
ful and accurate results. For a statistically stationary process, the time 
average of a quantity is independent of the placement of the interval over 
which the computation is made. This is expressed mathematically as 

4 b {J T f < > dt > ■ 0 < 9 > 

8 
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An equation for the fluctuating velocity components is obtained by 
subtracting Eq. 5 from Eq. 3. This yields 

XF ' (t) = a i U' (t) + b 1 V'(t) + c 1 W'(t) ' (10) 

The auto correlation of the fluctuating Doppler frequency F. ' (t) is 

T 

— f F ■ (t) F ,'(t+r)dt = X 2 F( ' (t) F. ' (tn) = 

J J 1 I J 1 

a. 2 U'(t) U ' (t+T ) + b:j 2 V'( t) V : (t+r) + c. 2 W'(t) W'(t+T) 

+ a^b^ U'(t) V 1 (t+x ) + a^c^ w'(t) u'(t+r) 

+ a.b. V'(t) U 1 (t+x ) + b.c. V'(t) W'U+t) (11) 

+ a.c. U'(t) W 1 ( t+T ) + b iC . W' (t) V' (t+x) 

For x = 0, Eq. 11 becomes 

X 2 F i , (t)2= ai 2 U'(t) 2 + b. 2 M' {t) 2 + c. 2 W'(t) 2 

+ 2a j bj U'(t) V'(t) + 2 a.c, U'(t') W'(t) 

1 1 1 1 (12) 

+ 2b iCi V'(t) W'(t) 


Equation 12 contains six unknowns, these being the variances and the 
mean instanstaneous cross products of the three velocity fluctuations. Thus 
in order to solve for the six unknowns, six independent measurements at six 

different geometrical configurations of the single LDV system would be required. 

Equation 11 provides a relation from which the auto and cross correlation 
functions may be obtained. This equation contains nine unknowns for each time 
delay. Thus, for stationary data, all the correlations may be solved using 
nine independent equations obtained by conducting LDV measurements at nine 
independent geometrical configurations. 
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It can further be shown that the spectra of. the individual time functions 
are related to the spectra of the correlation functions. Consider two functions 
and their respective Fourier representations 

oo 


f(t) = -2- 


g(t) = 


E 

n=l 

co 

E 


a cos nut +b sin n w t 
. n n 


c n cos n id t + d p sin n to t 


( 13 ) 


( 14 ) 


Where: 


Ti 

2_ 

2_ 


f 


Ji 


f(t) cos n w t dt 


Ji 


J f(t) sin n a) t dt 

O T 

J g(t) cos n w t dt 


A 


g(t) sin n w t dt 


( 15 ) 

( 16 ) 

( 17 ) 

( 18 ) 


The cross correlation of f(t) and g(t) now becomes: 

, /■ 

f(t) g(t + t) = — / f(t) g(t + t) dt 


V 0 = 


a c c o 


E 

n=l 


a c . b d 
n n + n n 


a n d n - c n b n 


09 ) 


cos nwT + 


sin nwx 
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When f(t) and g(t) are identical functions, the correlation is known as the 
auto correlation, and a n = fc> n and b n '"= d . Equation (T9) becomes: 

CO 

R^x) = f(t) f(t+x) = ] — !1_ + — — cosnwT (20) 

ft 4 2 

n=l 


Equations (19) and (20) are the Fourier transformation of the cross and auto 
correlations respectively. The spectral amplitude of the cross correlation 
function is directly related to the spectra of the individual functions. The 
absolute value of the spectra of the cross correlation is known as the cross 
power spectra of the functions f(t) and g(t). Also the spectral amplitudes of 
the auto correlation are the power spectral amplitudes of the original time 
function. Thus the power spectra and cross power spectra may be calculated 
from the auto and cross correlation functions. 
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SIMULTANEOUS EQUATIONS EVALUATIONS 

Although the previous developments theoretically demonstrate that the 
single LDV system can provide a technique for obtaining both three dimensional 
mean and statistical flow information in a stationary flow field, the ex- 
perimental feasibility of using the techniques needs to be examined. The basis 
of the described techniques, hinge on the solution of sets of independent 
linear equations containing either three, six or nine unknowns. In the solution 
of such sets of equations, small errors in the LDV measurements at any given 
geometrical system arrangement, can result in exceedingly large errors in the 
solution to the unknown values in the set of equations. In order to assess the 
magnitude of such errors, calculations were conducted for the three, six, and 
nine equations cases for possible geometrical arrangements of a typical reference 
scatter type single LDV system. 

Figure 2 shows a schematic of the geometrical instrument arrangement of 
a single component reference-scatter type LDV instrument. Note that the refer- 
ence laser beam is along the Z axis. For this arrangement, the two variables 
which define an independent measurement are the angles 0 and a. The angle 0 
defines the location of the plane formed by the reference and scatter beams with 
respect to the Y axis. The angle a is the scattering angle and is formed by the 
intersection of the reference and scatter beams. 

In the calculations conducted for the 3, 6, and 9 equation cases, the 
LDV system was assumed to be at different scatter angles and the optical pack- 
age, which collects the scattered light, was assumed to be rotated to different 
values of 0. For example, in one of the calculations, a was assumed to be 12 
degrees, while the instrument was rotated to 0 values of 0°, 120° and 240°. 


7 



The equations of interest in the calculations, are: 
3 Equation Case for Mean Values 


XF. (t) = a n . U(t) + b i V (t) + c. W(t) (5) 

where i = 1 — * 3 

6 Equation Case for RMS and Instanstaneous Cross Products 

X 2 F i '(t) 2 = a i U ' ( t ) 2 + b i 2 V'(t) 2 + c.2 W'(t) 2 
+ 23^. U' (t)V' (t) + 2a-} Cj U'(t)W'(t) 

+ 2b.c. V ' (t )W' (t) (12) 

where i = 1 — *6 

9 Equation Case for Auto and Cross Correlations 

X 2 F i 1 (t) F. ' (t+t) = a. 2 U' (t)U' (t+x) + b i 2 ~V' (t)V' 

+ c -} 2 W 1 ( t ) W 1 ( t+x ) + a i b i U'(t)V'(t+T) 

+ a^c^ U'(t)W'(t+x) + a^b} V , (t)U'(t+T) 

+ b^c. V'(t)W'(t+T) + a^. W '( t ) U '( t+r ) 

+ b iCi W'(t)V'(t+T) (11) 

where i = 1 — ► 9 

For the laser system and arrangement shown in Fig. 2, the coefficients 
in the above equations are: 

a. . = -sin 0.J sin 

b. = cos 0.: sin a. 

i ii 

c. - = cos a. -1 

1 i 

Also, for the LDV system considered, it was assumed that a frequency 
tracker was used to track the Doppler frequency. This instrument provides 
a DC voltage output E^(t) ^ F.(t) and a AC voltage output (t) a F^ (t). 
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The constant of proportionality between the output voltage and input frequency 
may be written as 

where D has dimensions of velocity per volt. 

For the calculations made here, D was assumed to be either 3.90 (m/sec)/ 
volt or 6.58 (m /sec)/volt, which represented a typical range of values which 
can be achieved with the frequency trackers described in Ref. 1. 

In order to determine the sensitivity of a LDV geometrical arrangement 
for a given set of measurements, it is necessary to set-up the coefficients 
of the simultaneous equations which defines the configurations. For the three 
equation case, these equations may be written as 

Ej (t) = A 1 U(tj + Bl W) + Cj“w(t) 

Mtj = A 2 U(t) + B 2 W + C 2 W(t) (21) 

Ej'(t) = A 3 m + B 3 V(t) + c 3 W(t) 

where A = aj/D, Bj = bj/D,---etc. 

The sensitivity of the measurements of Ej, E 2 , and E g to a determination 
of the three velocity components is obtained by simultaneously solving Eqs. 21 
for the velocities in terms of the measured voltages. This can be written as 

uTo = L x E^tj + M i E 2 (t) + N x E 3 (t) 

W) = +M a Mt) + N 2 ~Mt) ( 22 ) 

W(tT= L 3 E 1^) + M 3~E 2 (t) +N 3 E 3 ( t ) 

In Eq. 22 the sensitivity of the voltage measurement to a velocity 

component, i.e. AVelocity / a V oltage is given by the coefficients L., M , and 

i i 

. These coefficients are obtained by inverting the matrix formed by the 
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coefficients in Eq. 21. 

By using typical values of A,., B., C_. in Eq. 21 for different geometrical 

arrangements of the LDV instrument, values of , M., and N. were obtained. 

11 i 

The geometrical values used in determining A^ , B. , and along with the values 

of the matrix inversion are given in Table i. This table gives the geometrical 
arrangement of the cases considered which gave the best values of L.. , , and 

and the geometrical arrangement which gave the worst values. The best case 
was defined to be the geometrical arrangement which gave the lowests values 
for the output coefficients L.. , , and N. , and the worst case was defined 

to be the geometrical arrangement which gave the highest values for the output 
coefficients. The significance of the two cases can best be seen by noting that 
in the best case the output coefficients are of the order 10 1 while the output 
coefficients are of the order from 10 2 to 10 s for the worst case. From a 
measurement standpoint this means that a voltage measurement of t 1 x 10~ 3 
volts in the best case, will correspond to a velocity measurement of t 1 x 10" 2 
m/ft and for the worst case to - 1 x 10° to - 1 x 10 2 m/sec. Thus, it can 
be seen that measurements conducted using the geometrical arrangement of the 
worst case will impose a severe restriction on the precision a voltage measure- 
ment must be made in order to obtain an accurate velocity measurement. 

After examining the results of the different cases that were considered 
for the 3 equation case the following trends were noted. 

1. Better accuracy is achieved when a is increased. 

2. Better accuracy is achieved when the three measurements are conducted 
at 120° intervals of 0. 

3. Better accuracy is achieved as the value of D is reduced. 

The six and nine equations which defined the different statistical 
quantities were also examined for different geometrical arrangements of a 
single component reference scatter LDV system. Using the equations given by 
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3. The auto and cross correlations for a three dimensional flow can be 
measured using nine independent measurements. 

The defining equations were examined for a single component forward 
scatter LDV system for different geometrical arrangements. It was found • 
for the arrangements considered; 

1. that the three mean velocities could be determined with a reasonable 
accuracy and 

2. that the higher moment statistical quantities would be very difficult 
to determine accurately. 

Only a limited number of geometrical arrangements of a particular LDV 
system were considered here. Any geometrical arrangement considered should 
always be guided by the physical constraints imposed by the LDV system and 
the flow facility, since many arrangements that are theoretically possible 
are not physically practical. 

It should be noted that reducing the flow field to a two dimensional 
case causes the required number of independent configurations to be noticeably 
reduced for the determination of any particular statistical quantity. That 
is only two configurations are needed for mean measurements, three for RMS and 
the cross product, and four independent observations to reduce all time cor- 
relations. It is of further interest to note that if the LDV signals are 
collected from three independent directions simultaneously (Ref. 2), the 
three demensional statistical quantities may be obtained directly with a 
minimum restriction on the precision with which the Doppler frequencies must 
be detected. 
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(10) and (11) the matrix of the coefficients were inverted to obtain the output 
coefficients which related the voltage measurement to the different statistical 
velocity, parameters. This investigation indicated that the output coefficients 
for the 6 equation case were of the order 10" 3 to 10 12 (m/sec)/vo1t and for the 
9 equation case of the order 10 1 to 10 X1 (m/sec)/vo1t. The magnitude of the 
output coefficients would require a extremely precise voltage measurement 
capability in order to obtain a satisfactory accuracy in the measurement of 
the statistical velocity parameters. Tables II and III provide a tabulation 
of geometrical arrangements considered for the 6 and 9 equation cases and the 
magnitudes of the output coefficients. 


CONCLUSIONS 

The general equations defining the three dimensional mean and fluctuating 
velocity components measured by a single LDV instrument were presented. It 
was shown that mathematically, a single component LDV system has the ability 
to measure both three dimensional mean and fluctuating velocity components 
assuming that the flow field is stationary. In particular, it was found that: 

1. The three component mean velocities can be measured by conducting 
LDV measurements at three independent geometrical configurations 
of a LDV system. 

2. The RMS and instantaneous cross product averages for a three dimen- 
sional flow can be measured by making LDV measurements at six inde- 
pendent geometrical configuration. 
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Table 1 

Results, of Sol utlon for 
Three Independent LDV Equations 

Best Case Considered: D = 12.8(ft/sec/vol t). 

otj = a 2 = a 3 = 28° 

0 j = 0°, 02 = 120°, 0 3 = 240° 

Input Coefficients (volts/(ft/sec) 

Aj = 0.0 = 3.643 x 10' 2 C = -9.084 x 10' 3 

A 2 = -3.155 x 10" 2 B 2 = -1.821 x 10 -2 C = -9.084 x 10~ 3 

A 3 = 3.155 x 10- 2 B 3 = -1.821 x 10' 2 C 3 = -9.084 x 10' 3 

Output Coefficients ((ft/sec)/volt) 

Li = -5.769 x 10- 4 Mj = -1.584 x 10 1 Nj = 1.584 x 10 1 

L 2 = 1.829 x 10 1 M 2 =-9.148x10° N 2 =-9.158x10° 

L 3 = -3.669 x 10 1 M 3 = -3.669 x lO 1 N 3 = -3.669 x 10 1 

Worst Case Considered: D = 21 .6( (ft/sec)/volt) 

ql 1 = a 2 = a 3 = 5° 
e 3 = 270°, 0 2 = 280°, 0 3 = 290° 

Input Coefficients (volts/(ft/sec)) 

A, = 4.027 x 10 -3 B 2 =2.165xlO- B Cj = -1.758 xlO' 4 

A 2 = 3.966 x 10 -3 B 2 = 6.993 x TO" 4 C 2 = -1.758 xlO' 4 

A 3 = 3.784 x 10 -3 B 3 =1.377xl0- 3 C 3 = -1.758 xlO' 4 

Output Coefficients ((ft/sec)/volt) 

L x = -7.923 x 10 3 •* 1.609 x 10- . = -8.172 x 103 

L 2 = -2.123 x 103 ‘M 2 = 2.838 x 103 N 2 = -7.150 x 10* 

L = -1.871 x 10 s M = 3.686 x 10* N, = -1.871 x 10* 

3 3 J 
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Table. II 


Geometrical Arrangements Considered 
for Solution to Six LDV Equations Defining the 
Three Dimensional RMS and Cross Product Terms. 


a* 

0! 

0 2 

0 3 

04 

e 5 

9 6 

Output Coeff. 

Magnitude 

(ft/sec)/volt)) 

5 

0 

45 

90 

135 

270 

315 

O 

1 

i-* 

1 

o 

M 

IO 

8 

0 

45 

90 

135 

270 

315 

10-M0 11 

9 

0 

45 

90 

135 

270 

315 

10" 2 -10 12 

10 

0 

45 

90 

135 

270 

315 

04 

O 

1 

OJ 

O 

11 

0 

45 

90 

135 

270 

315 

10 3 -10 16 

12 

0 

45 

90 

135 

270 

315 

0-10 16 

5 

0 

60 

120 

180 

240 

300 

10- 3 -10 12 

12 

0 

60 

120 

180 

240 

300 

lO-MO 11 

5 

270 

280 

290 

300 

310 

320 

10 6 -10 13 

12 

270 

280 

290 

300 

310 

320 

r-* 

o 

1 

J- 

o 

5 

270 

285 

300 

315 

330 

345 

10 s -10 12 

8 

270 

285 

300 

315 

330 

. 345 

*■4 

o 

1 

o 

11 

270 

285 

300 

315 

330 

345 

10 5 -10 12 

12 

270 

285 

300 

315 

330 

345 

10--10 11 


* All angles in degrees 
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Table I LI .. 

Geometrical Arrangements Considered for 
Solution to the Nine. LDV. Equations Defining the 
Three Dimensional Auto and Cross Correlations 


a* 

0 i 

0 2 

6 S 

0 4 

0 5 

0 6 

6 7 

9 e 

0 9 

Output Coeff. 

Magnitude 

(ft/sec)/volt)) 

5 

0 

30 

60 

90 

120 

150 

240 

270 

300 

10 3 -10 12 

8 

0 

30 

60 

90 

120 

150 

240 

270 

300 

10 3 -10 n 

9 

0 

30 

60 

90 

120 

150 

240 

270 

300 

0 

u> 

1 

O 

*-» 

10 

0 

30 

60 

90 

120 

150 

240 

270 

300 

10 4 -10 n 

11 

0 

30 

60 

90 

120 

150 

240 

270 

300. 

10 4 -10 1J 

12 

0 

30 

6 

90 

120 

150 

240 

270 

300 

10 4 -10 12 

5 

0 

40 

80 

120 

160 

200 

240 

280 

360 

10 3 -10 12 

12 

0 

40 

80 

120 

160 

200 

240 

280 

360 

10 3 -10 12 

5 

270 

280 

290 

300 

310 

320 

330 

340 

350 

10 4 -10 12 

8 

270 

280 

290 

300 

310 

320 

330 

340 

350 

1 0 4 -l o 1 1 

12 

270 

280 

290 

300 

310 

320 

330 

340 

350 

10M0 11 

45 

270 

280 

290 

300 

310 

320 

330 

340 

350 

0 

fo 

1 

o 

CD 

90 

315 

325 

335 

345 

355 

5 

15 

25 

35 

10 3 -10 7 


* All angles in degrees 
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Dual Beam Forward Scatter System 

X 
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FIGURE 2. SCHEMATIC OF A TYPICAL ANGULAR ARRANGEMENT OF A 
SINGLE-REFERENCE-SCATTER LDV SYSTEM 
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